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Abstract. We note that a tridiagonal matrix representation of the algebra of the 
partially asymmetric exclusion process (PASEP) lends itself to interpretation as the 
transfer matrix for weighted Motzkin lattice paths. A continued fraction ( " J-Fraction" ) 
representation of the lattice path generating function is particularly well suited to 
discussing the PASEP, for which the paths have height dependent weights. We show 
that this not only allows a succinct derivation of the normalisation and correlation 
lengths of the PASEP, but also reveals how finite-dimensional representations of the 
PASEP algebra, valid only along special lines in the phase diagram, relate to the 
general solution that requires an infinite-dimensional representation. 
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1. Introduction 

Although the asymmetric exclusion process (ASEP) — a model in which hard-core 
particles hop in a preferred direction along a one-dimensional lattice — has been 
reinvented in various different guises over the years, it is only relatively recently that 
exact solutions for the steady state (s) of the model have been available. The solution 
of the ASEP with open boundary conditions in [1] using a matrix product ansatz was a 
landmark in the study of driven diffusive systems. 

As discussed in a recent review of the matrix product approach to solving for 
the steady state of nonequilibrium Markov processes [2], there are a range of different 
methods for analysing the thermodynamic phase behaviour of the simplest versions 
of the ASEP. By contrast, more general models — collectively known as the partially 
asymmetric exclusion process (PASEP) — that admit particles to hop in both directions 
in the bulk, and even more generally to enter and exit at both left and right boundaries, 
have so far been studied only through a diagonalisation of the matrices appearing in 
the formalism [3-5]. In this work, we extend a technique that previously admitted 
an extremely quick derivation of the ASEP phase behaviour under various updating 
schemes [6,7] to these more general models. 

The idea is to consider the behaviour of a "grand-canonical partition function" for 
the model. More precisely, we examine the generating function of the normalization of 
the nonequilibrium steady-state distribution over an ensemble of different lattice lengths 
whose mean is controlled by a fugacity. The thermodynamic phase behaviour can then 
be read off from the singularities of this generating function. Whilst obtaining this 
generating function is straightforward for the ASEP [6,7], a convenient closed form for 
the PASEP has remained elusive. 

Our aim here is to demonstrate that a representation of the generating function that 
allows the thermodynamic phase behaviour to be determined with relative ease takes 
the form of an infinite continued fraction. This we arrive at through an interpretation of 
the PASEP normalization as the (equilibrium) partition function of lattice paths, which 
we discuss in Section [3] after recalling the model definition and its basic properties. In 
Section [4] we show how to analyse the singularities embedded in the continued fraction 
representation. The results we obtain are, of course, equivalent to those obtained within 
other approaches [3-5]. However, given that continued fractions are not frequently 
encountered in statistical mechanical contexts, we feel there is some value in using the 
PASEP as an illustrative example of how to handle them. 

We find that the analysis is intimately related to an approach based on finite- 
dimensional matrix representations [8,9], exact along special lines in the phase diagram, 
and that the continued fraction shows how these particular solutions and the general 
solution are related. We further show that the continued-fraction approach extends 
to the most general version of the PASEP, solved in [5], and that one can access both 
currents and correlation lengths through it. Finally, we return to the lattice path picture 
to elucidate the equilibrium counterpart of a nonequilibrium phase transition identified 
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in [4] that occurs when the bias on bulk hop rates opposes that imposed by the boundary 
conditions. 

2. Model Definition and Basic Properties 

The dynamics of the PASEP take place on a finite one-dimensional lattice with open 
boundaries. In its simplest form, the microscopic dynamics of the PASEP are specified 
by four rates, one of which can be set to unity by an overall scaling. For a rate 
A associated with a particular event, the probability that the event happens in an 
infinitesimal time interval At is XAt. Moves that would lead to two particles occupying 
a single lattice site at any one time are prohibited due to the hard-core repulsion between 
them. 

In the PASEP particles are inserted onto the left boundary site (when empty) at 
a rate a and removed from the right boundary site at a rate f3 , see figure [T} Once on 
the lattice a particle hops by one site to the right at rate 1 or by one site to the left at 
a rate q when sites are available (i.e. empty). It is possible to expand this set of moves 
to allow particles to enter at the right at a rate S and exit at the left at a rate 7, while 
still retaining the solvability of the model [5]. 

In all the models we consider we want to calculate Z, which normalises the statistical 
weight, f{C), of a lattice configuration, C, in the steady state. This is given by 

Z = J2f(C)i (2-1) 

c 

so the normalized probability of being in state C is P(C) = f(C)/Z. 

The weights themselves are obtained through the stationarity condition on the 
transition rates W(C — > C), 

J2 [f(C')W(C' - C) - f(C)W(C -> C')] = , (2.2) 

where W(C — > C) is the probability of making the transition from configuration C to 
C in a single timestep. This is less restrictive than the detailed balance condition for 



equilibrium states, which is obtained when the sum in equation (2.2) vanishes term by 
term. 

The solution of the ASEP in [1] and the PASEP in [3-5] made use of a matrix 
product ansatz [2]. In this the steady-state probability P{C) of a configuration of 
particles C on a chain of length N is represented by an ordered product of matrices 
X\ X 2 . . . X N where X^ — D if site i is occupied and X, = E if it is empty. We expect 

A A a JL A JL A 



Figure 1. Typical particle configuration and allowed moves in the PASEP model. 
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P(C) to be a function of both the number and position of particles on the lattice, which 
suggests the choice of non-commuting objects, matrices, for the ansatz. To obtain a 
scalar probability value from this matrix product it is sandwiched between two vectors 
(W\ and \V): 

P(C) = . (2.3) 

z N 

The factor Z^ is included to ensure that P(C) is properly normalised. This latter 
quantity plays the role of a partition function in equilibrium problems 

Z N = (W\(D + E) N \V) = (W\C N \V) , (2.4) 

where we have defined C = D + E. Indeed, we shall see in what follows that is the 
partition function for an equivalent two-dimensional lattice path problem. 

The algebraic properties of the matrices D and E can be deduced from the master 
equation for the dynamics of the ASEP, PASEP and various other related models [2]. 
For the variant of the PASEP discussed above, sufficient conditions for equation (2.3) 
to hold are 

DE — qED = D + E , (2.5) 
a(W\E = (W\ , (2.6) 
/3D\V) = \V) . (2.7) 

These relations allow one to calculate Z^ and other quantities of physical interest by a 
range of methods, such as "normal-ordering" of the matrices, or through use of explicit 
representations [2]. 

In this work, we focus on an approach based around the generating function 
of Zn, namely Z{z) = Y2n %nz n , which can be thought of as a "grand-canonical" 
normalization. As is well known [10], the large- N form of the "canonical" normalization 
(Zn) can be determined from the dominant singularity z cr of Z{z). Typically, Zn ~ 
z~ r N N~ u where the exponent v > depends on the nature of the singularity. Then, by 
defining a "reduced free energy" / via 

/ = - lim ^lnZ N , (2.8) 

we find / = lnz cr . Nonanalyticities in / can then be associated with phase transitions 
in the physical system [2, 6, 7, 11]. 

For orientation, let us recall the results for the ASEP, which has q = 0. The 
canonical normalization can be shown by direct matrix reordering [1] to be 

7 _^ p(2N-i- P )\(i/^-(i/ a y^ 

N ^ M(N-p)\ (l/P)-(l/a) ■ { ■ ] 

Performing the summation [6] gives the grand canonical normalization 

Z{Z) = (a - x(z))(/3 - x(z)) ' (2,10) 
where x(z) = (1 - VI - 4z)/2 . 



Continued Fractions and the Partially Asymmetric Exclusion Process 



* 

p 


i 






LD 


MC 


/ 






2 








/ 


HD 



a 



Figure 2. The phase diagram of the ASEP. Here HD, LD and MC denote the 
high-density, low-density and maximal-current phases, respectively. 



This function has a pole at x(z) = a when a < \ and similarly at x(z) = (3 when 
(3 < \. These correspond to z cr = a(l — a) and z cr = (3{1 — (3). When neither of these 
poles contribute, all that remains is the square-root singularity at z cr = |. This allows 
one to very quickly establish the behaviour of the reduced free energy as a function of 
a and (3: 



f 



(2.11) 



In [1] for a, (3 > 1/2 

ln[a(l-a)] for (3 > a, a < 1/2 . 
ln[/3(l -p)] for a > (3, f3 < 1/2 

It turns out that for the ASEP, z cr corresponds to the particle current, and x(z cr ) 
the bulk density in the thermodynamic limit. Hence, the phase diagram for the 
model, Figure [2j is quickly recovered using this generating- function (or grand-canonical) 
analysis. Further details of these methods as applied to the ASEP can be found in [2,6,7]. 
In the remainder of this work, we show how to elicit the structure of the grand-canonical 
normalization of the PASEP, where direct summation of the canonical normalization, 



given explicitly in [4], does not lead to a compact expression like (2.10). 



3. Explicit Matrix Representation and Lattice Path Interpretation 

A useful route to the grand-canonical normalization for the PASEP is via the generating 
function for an ensemble of lattice paths, which in turn can be read off from an explicit 
representation of the matrices and vectors appearing in the equations (2.5), ( |2.6 ) and 



(2.7) that define the matrix algebra. A number of representations are known, see [1-3,5]; 
the one that is of use here is that for which the vectors (W\ and | V) have nonzero entries 
only in their first element: 



( Wq \ = (W\ = 0, 0, • ■ ■) \V q ) = \V) = h^\l, 0, 0, • • -Y , (3.12) 

where h is a constant to be given shortly. One can verify that, with this choice of 



boundary vectors, the following tridiagonal representations of D and E satisfy (2.5) 
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dZ6b and dZ7h: 



,0, 



1-9 



1-q 



( 1 + P 

1 + /3g ^ 

1 + pq 2 

( 1 + 5 

■y/ci 1 + 5g 

v /c^ 1 + aq 2 

X '■■ 



... \ 



(3.13) 



The various parameters that appear are 
1-q 



a 



P 



- 1 



a 
1-q 

P 



c n = (l-q n )(l-a(3q n - 1 ), 
1 A (ap) n 



h () 



(3.14) 

(3.15) 
(3.16) 
(3.17) 



in which we have used the standard notation for (shifted) g-factorials 

n-l 
3=0 

(a;g)o =i, 

(a, 6, ... , c; g) n = (a; g) n (6; g) n . . . (c; g) n . (3.18) 

There are various ways to arrive at an interpretation in terms of lattice paths from 
a matrix representation. One was suggested by Brak and Essam [12, 13], who used the 
fact that D + E = DE for the ASEP to interpret the D and E as odd-even and even- 
odd height transfer matrices separately. In later works [6,7], the path interpretation 
was inferred from the grand-canonical normalization once this had been obtained by 
another means. Here, we shall take the most direct approach, which is to associate a 
height n > above the origin with the vector \n) = (0 • • • 1 • • -) T (i.e., n is the 
number of zero entries that appear before the single nonzero entry). We then interpret 
(m\X\n), where X is some combination of D and E matrices, as the weight of paths 
connecting a point at height n to another point at height m. 

Of particular importance is the matrix X = C , which appears in Eq. (2.4) for the 
normalization. From the above expressions we have that 

/ 2 + ti + p ^c[ ••• \ 

I ^fc~\ 2 + (a + P)q ^[02 

1-q ^/ci 2 + (a + (3)q 2 '•• 



C = D +E a 



V 



.(3.19) 
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h + + + +^+ + + + H h + + + ^ ^ H 

Figure 3. Two Motzkin path transliterations of the tridiagonal matrix representation 
of the PASEP matrices. In both cases, up-step down-step pairs contribute a weight c ni 
where n is the height above the origin of the upper end of the steps. In the left-hand 
figure, there is a single type of horizontal step that contributes a weight d n , where n is 
the height of the segment above the origin. The right-hand figure shows an equivalent 
interpretation, in which there are two types of horizontal steps, shown with diamonds 
and circles at the their midpoints. These contribute weights and d„, respectively. 



The matrix element (m\C N \n) then gives the combined weight of paths that begin at 
height m, end at height n and contain N steps, each of which may raise or lower the 
height by one unit, or leave the height unchanged. Since the height is a nonnegative 
quantity, n > 0, these paths may never descend beneath the origin. Additionally, since 
(W\ oc (0| and \V) oc |0), the paths that contribute to Zn begin and end at the origin. 
Paths with these properties are known as Motzkin paths. 

The weight of various path components can now be obtained by inspecting the form 
of C q . For the path to begin and end at the origin, every up-step must be accompanied 
by a down-step; each up-step down-step pair connecting height n — 1 to height n 
contributes a weight c n to the path. Each horizontal step at height n contributes a 
weight d n = 2 + (5 + (3)q n to the path. An alternative interpretion has two types (or 
"colours" ) of horizontal path segments, one of which contributes a weight = 1 + aq n 
and the other d b n — 1 + (3q n . See figure [3} To arrive at the canonical normalization 
Z/v, we sum over all paths of length N, and multiply by the factor h /(l — q) . The 
grand-canonical normalization is constructed by summing paths of all lengths, weighting 
each segment by z, and finally multiplying by h Q . In the next section, we shall see an 
equivalent recursive construction which can be expressed as a continued fraction. 

As a check of the path representation one can compare the expression for which 
emerges from directly evaluating the matrix product expression for the stationary state, 

Z N = (W\V)( — \ J2 R NM B n&fc<l), (3-20) 

\ ^/ n =0 

where 

rnm = zlTh-ir i( N z 2k ) - (vJ-V 2 )] ^ m , 

(3.21) 
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and 



n 

B n (a,P; q) = [ 



k=0 



n 



j3 



(3.22) 



with the weights which emerge from the contributing paths at low orders. In the above 
we have used the standard notation for the g-binomial coefficient 

(q;q) n 



Ikl 



(3.23) 



(q;q)n-k(q;q)k ' 

For instance, taking the simple case of Z 2 a direct calculation using the above formulae 
gives 



5 - q + a 2 + a/3 + (3 2 + qa(3 + 4a + 4(3 



(3.24) 



(i- q y 

where we have dropped the overall normalization (W|V). The (1 — q) 2 denominator 
disappears when this is written in terms of the original rate parameters a and (3, 

a(3q + a 2 + (3 2 a + a(3 + a 2 (3 + (3 2 



a 2 f3 2 



(3.25) 



The expression for Z 2 in (3.24) can be seen to be the sum of weights for an 
up/down step pair from level to 1, given by (1 — q)(l — a(3)/(l — q) 2 and all four 
possible combinations of two horizontal steps at level zero of either "colour" , given by 
(2 + a + (3) (2 + a + 13) /{I — q) 2 (where we have again dropped the overall normalization 
(VK|V)). While the diagrammatics becomes increasingly complicated for larger Z^, the 
principle remains the same. 



4. The Continued Fraction Representation of the Path Generating Function 



Despite the availability of an exact expression for the canonical normalization, Eq. (3.20 ) 
from [4], we have not been able to find a convenient expression for its grand-canonical 
counterpart due to the g-dependence of the PASEP weights. A more fruitful route is 
to represent the lattice-path generating function as a continued fraction, a procedure 
first expounded by Flajolet [14]. As we now show, this representation can be read off 



more-or-less directly from the matrix representation (3.19), and the form that emerges 
is particularly well adapted to discussion of the PASEP. 



First, for convenience, we subsume the prefactor 1/(1 — q) appearing in (3.19) into 
the parameters 

2 + (6 + (3) q n 



d n 

L-q 

h _ n n \(^ - ™Rn n - l \ 

(4.26) 



l-q 
(l-q) 2 
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so that then 



C„ 



d 




\JC\ 





d 2 



\ 



\ 



J 



(4.27) 



In the lattice-path language, this means that horizontal steps at height n are weighted 
by d n and up-down step pairs between heights n and n + 1 are weighted by c n . 

Let now Ai n (z) be the generating function of weighted Motzkin paths that start 
and end at height n, never go below this height, and have their lengths N > counted 
by powers of z. That is, the coefficient of z N in A4 n (z) is the weight of such paths 
of length N. The grand-canonical normalization for the PASEP is then given by 
Z{z) = hoM-o(z). Let us suppose that M. n +i(z) is known for some n > 0. Then, 
we can construct M. n (z) by concatenating contiguous components of two types: (i) 
sequences of horizontal segments of arbitrary (possibly zero) length at height n; and (ii) 
Motzkin paths starting at height n + 1 enclosed by an up-down pair. Denoting these 
components schematically as and /M. n +i \ respectively, we can write the recursion 



M n = _ + _/M n+1 \_ + _/M n+1 \^/M n+ i\_ +■■■ (4.28) 
= _ (1 + [/M n+l \_] + [/M n+1 \_] 2 + ■■■) (4.29) 

The generating function for a (possibly empty) sequence of horizontal segments, each 
weighted by d n , is simply (1 — d n z)~ x . An up-down pair / • • • \ from height n to n + 1 
contributes the weight c n+ iz 2 . We thus arrive at the generating- function recursion 

(i-LzY 1 i 



M n (z) = 



.(4.31) 



1 - c n z 2 M n +i(z)(l - d n z) 1 1 - d n z - c n+1 z 2 M n+ i(z) 

Starting at n = and iterating, we find that the generating function Z(a, (3, q, z) 
for Motzkin paths of arbitrary length (and hence the grand-canonical PASEP 
normalization) is given by the infinite continued fraction 

1 



Z(a,(3,q,z) = 



1 — d z 



C\Z 



(4.32) 



1 — d\z 



c 2 z 2 



1 — d 2 z — 



c 3 z 2 



where we have dropped the factor h since this does not contribute to any physical 
quantities. Such a continued fraction containing both z and z 2 terms is usually denoted 
a Jacobi continued fraction, or "J-fraction" for short [15]. 

Before considering the case of general q let us take q = and see how the expression 
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for the grand canonical normalization, Eq. (2.10), is recovered. We have 



Z{aj,0,z) = , (4.33) 

1 — d()Z — 2 

l-2z - 



1 - 2z 



z 2 



1 - 2z 



z 2 



where d = l/a + 1//3 and C\ = 1/a + — l/(a/3) = k 2 when q = 0. Note that in 
this case, the continued fraction is periodic after the first level. That is, in the above 
notation, 

Hence, we must have that A4i(z) = M%{z) = ■ ■ ■ and hence 

Mi(z) [I -2z-z 2 Mi(z)] = 1 (4.35) 

or 

l-2z- VI -4z 
~2? 

which is the generating function familiar from many Catalan counting problems. For 
n = 0, we have 



Mi(z) = ^ , (4.36) 



Z(a, 0, 0, z) = M (z) = — \ , (4.37) 

1 — d z — c\Z z M.\\z) 



which coincides with (2.10) when both are expanded and rationalized. 



At q = the luxury of being able to sum the continued fraction to get equation 



(4.37) makes the phase structure, which (as previously discussed) is determined by the 
singularities of Z(&, /3, 0, z) in z, immediately apparent. If we are not able to easily sum 
explicitly the continued fraction, as is the case for the PASEP, we can use more indirect 
methods to determine the singularities. We focus here on the "forward-bias" regime, 
q < 1, in which the continued fraction has a finite radius of convergence in the complex- 
z plane. In the reverse-bias regime, q > 1, the continued fraction is unconditionally 
divergent, a fact we will interpret physically in Section [7j 

We first appeal to Worpitzsky's theorem on the convergence of continued fractions 
[15, 16] which states that a continued fraction of the form 

l -a 2 (4-38) 



1 + 



a 3 



04 



converges if the partial numerators a p satisfy 

\a p \ < 1/4, p = 2,3,4, ... . (4.39) 



Continued Fractions and the Partially Asymmetric Exclusion Process 



11 



For Z(z) given by (4.33) this translates to a radius of convergence z cr given by 

= 1 Vn 



(1 - rf n _i2 cr )(l - 
and shows that 

z cr -»• (1 - g)/4 



dnZcr) 



(4.40) 



(4.41) 



as n — > oo. To decide if this is the dominant singularity, one must also divine the 
location of any poles in the complex- z plane from the continued fraction (4.32). 

The strategy is to examine the n th convergent of the continued fraction, that is, 
the expression obtained by truncating the continued fraction at the n th level (counting 
from zero). Denoting this as K n , we have that 

1 



Kn 



1 — dnz 



1 — dnZ — 



1— d\z 



(4.42) 
(4.43) 



and so on. We observe that the continued fraction (4.33) is given exactly by the 
convergent K n if c n+ i = 0. Inspection of (4.26) reveals that this occurs on the special 
line in the phase diagram given by 

a(3 = q~ n . (4.44) 

We observe that on such special lines, the matrix C q decomposes into two blocks: the 
first is (n + l)-dimensional, whilst the second does not contribute to the normalisation 
because only the first elements of the boundary vectors are nonzero. Finite-dimensional 
representations of the matrix algebra for the PASEP were used prior to the advent of the 
full solution to study the model along these lines, and to conjecture the phase behaviour 
elsewhere [8,9]. 

To locate the poles of K n in the complex-2; plane, we use the fact [17] that the 
reciprocal of each convergent K n can be expressed as 



B r . 



(4.45) 



Bn 



(4.46) 
(4.47) 

1. Hence, 



where A n and B n both satisfy the same three-term recurrence 

A n (z) = (1 - d n z)A n „ 1 (z) - c n z 2 A n _ 2 (z) , 
B n (z) = (1 - d n z)B n ^ l (z) - c n z 2 B n _ 2 {z) , 

but with slightly different initial conditions: A_ 2 = B-\ = and A_i 
K n has poles that coincide with the zeros of A n (z). These relations for A n (z) and B n {z) 
hold for arbitrary choices of at and (3. Let us now define the function A n (z) that is 
given by A n (z) for the particular choice (3 = l/(q n a). That is, the zeros of A n {z) give 
the poles of the PASEP grand-canonical generating function along the line in the phase 
diagram along which the (n + l)-dimensional matrix representation is exact. 
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As a consequence of the recurrence (4.46) the A n (z) are closely related to the al- 



Salam-Chihara polynomials [18] and may be written explicitly as 

q) k (ae~ ld ;q) k 



A n (z) 



ae ie - 



[1 - q) n + l a n+l 



k=0 



a(3;q) k (q;q) k 



(4.48) 



where cos(#) = (1 — q)/ (2z) — 1. If we take (3 = l/(q n a) to obtain A n (z) all the terms in 
the sum in (4.48) will be zero apart from the (n + l)th since they contain a (1 — a(3q n ) 



factor in the numerator, which is cancelled only in the last term. Using 



(ae ld - 



j=0 



we thus find 



where 



(l + aq k 



aq" 



(4.49) 



(4.50) 



(4.51) 



Equally, we could have taken a = l/(q n /3) and obtained the same result with a «-> j3. 
These roots correspond to the reciprocal of the eigenvalues of the finite-dimensional 
matrices C q examined in [8, 9] (note, however, that the representation used in those 
works differs from that used here). This fact can be understood by making the ansatz 



z n \a 



(n) 



where Cq^ is the submatrix formed by the first n rows and 



A n (z) 

columns of C q . After substituting into (4.46), one obtains the expression that results 
from a co- factor expansion along the bottom row in the determinant |C™ — A|, as long 
as A is identified with -. 

z 

We see, then, that the sequence of generating functions A n (z) has a truly remarkable 
property: all poles of A n {z) are poles of A m (z) where m > n. In principle, a completely 
different set of singularities could have been obtained after truncating at a deeper level. 
Furthermore, when q < 1, \z n \ < \z m \ if n < m, and hence all convergents have the same 
dominant pole, zq. Writing this in terms of the original rates a and (3 gives 

1 — q a(l — q — a) 



a 



l/a 



or 



l-q 

13(1 -q- P) 



(4.52) 



(4.53) 



2 + (3 + 1/(3 l-q 

Taking now the number of levels (or, equivalently, the dimensionality of the matrix 
representation) to infinity, we find that the dominant singularity will always be the 
smallest of (4.41), (4.52) or (4.53). This observation allows us to recover the reduced 
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free energies for the PASEP previously calculated in [3,4], and identify the regions of 
the phase diagram within which they apply: 

'1 



/ 
/ 
/ 



In 



In 



In 



Q 



a(l — q — a) 



0(1 



- q 

Q-0) 



1-q 



for a, (5 > (1 - q)/2 , 
, for (3 > a, a < 
, for a > (3, (3 < 



1 



ff)/2, 
?)/2- 



(4.54) 



As for the ASEP, this phase structure comprises high- and low-density phase and a 
maximal current phase. The key differences lie in the g-dependent position of the 
second-order transition line, a = /3 = (1 — q)/2, and in the subextensive corrections to 
the free energy, which manifest themselves in the density profile at finite distances from 
the boundaries [20]. 

As noted above, this analysis of the PASEP phase behaviour is similar to that 
based on finite-dimensional matrix representations [8,9]. The main additional benefit 
of the present approach is that the relationship between the general expression for the 
normalisation, given through the non-terminating continued fraction (4.32), and the 
versions that apply only on special lines is clearer. For example, the singularity that 
arises from non-convergence of the continued fraction, and which governs the maximal- 
current behaviour, is not detectable by inspection of the finite dimensional matrices. 
Furthermore, some light is shed on the "continuity arguments" expounded in [8,9] that 
lead one to extend physical properties found along special lines to regions of the phase 
diagram, in that each additional level of the continued fraction yields a singularity that 
is subdominant to all of those present at the previous level. We remark that the fact 
that the asymptotically dominant singularity is present even in the one-dimensional 
representation is suggestive of an explanation as to why the phase diagram and currents 
obtained from mean-field theory coincide with the exact results. 



5. Lattice Paths for Non-Zero "Wrong Direction" Removal and Injection 
Rates 7 and 5 

As we noted in the introduction it is possible to generalize the PASEP boundary 
conditions to allow "wrong direction" injection and removal rates. In addition to 
injection at a rate a at the left boundary particles may also be removed there at a 
rate 7. Similarly, in addition to removal at a rate f3 at the right boundary particles may 
also be injected there at a rate 5. 

A tridiagonal representation of D and E with these generalized boundary conditions 
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still exists for the case of non-zero 7 and S [5] and which can be written as 



/ 1 + d\ d[ 



D n = 



1-5 



\ 



4 



1 + d\ d\ 
d\ 1 + d\ 



E n = 



1-q 



<wg = /4 /2 (i,o,o,---), |y,) = ^(1,0,0,...) 



/ 1 + eJ 



V 



l + e\ 



1 + e 



/ 



1/2/ 



(5.55) 



(5.56) 



where the expressions for $ n , d^, e^, e h n (a notation borrowed from [5]) and h are given 
in the appendix. 

We can define C q = D q + E q in an analogous manner to the earlier discussion which 
leads to a generating function of the form 

Z(a, b, c, d, q, z) = 



1 — doz 



C\Z 



(5.57) 



1 — d\z 



c 2 z* 



1 — doz — 



C3Z 



where the coefficients are now 

d„ 



2 + 4 + 4 



1-q ' 

. (4 + 4K4±^) (* ^ 

C n +l ~ ^ _ ^ 2 • 

Even with these more complicated coefficients the singularity structure of the generating 
function is identical to the 7 = 5 = case when q < 1. We can still extract poles in 
order to discern the high- and low-density phases. Explicitly, 



~ (l-q n - 1 abcd)(l-q n + 1 )(l-q n ab)(l-q n ac)(l-q n ad)(l-q n bc)(l~q n bd)(l-q n cd) 
° n + 1 ~ (l-q 2n - 1 abcd)(l-q 2n abcd) 2 (l-q 2n + 1 abcd)(l-q) 2 

(5.59) 

(in which the parameters a, b, c and d are also given in the appendix) so that we can 
pick up the various poles depending on the zeros of the numerator. For non-zero 7 and 
5 Askey- Wilson polynomials [19] play the role of the al-Salam-Chihara polynomials so 
in this case the denominators of the convergents are given by 



(ab,ac,ad;q)n+iz n+1 ;q) fc (q"afecd)(q e^;q) k (a e -^;q) k fc 

(l- g )n+l a «+l Z^fc=0 (ab;q) k (ac;q) k (ad;q) k (q;q) k « 

(5.60) 
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where cos(#) = (1 — q)/(2z) — 1 as before. If we now choose q n ab = 1, for example, to 
enforce c n+ \ = we find that only the final term in the sum contributes again, resulting 
in the dominant poles 

l-q 



2 + a + l/a 



(5.61) 



or 

Z " = 2 + b + 1/6 (5 - 62) 
depending on the relative size of a and b. The maximal current phase is still observed in 
a similar manner to the 7 = 5 = case via Worpitzky's theorem but a major difference 
is that the reverse bias phase is absent when q > 1, since the particles can now escape 
at both ends. Algebraically this is manifested in the reflection symmetry relating the 
parameters for q < 1 and q > 1 when 7 and 5 are non-zero: 

a, b, c,d,q — > b~ x , a -1 , gT 1 , c _1 , . (5.63) 

We remark that this continued-fraction representation provides a very quick route to 
the identification of the dominant singularities in the generating function, and therewith 
the extensive part of the free energy. 



6. Correlation Lengths 

We have so far restricted our discussion to the normalization itself but it is also 
possible to extract the correlation lengths, which determine the various sub-phases 
within the PASEP phase diagram, in a rather straightforward manner. We have seen 
that the leading singularities of the grand-canonical normalization Z determine the 
phase structure of the PASEP and that these can be extracted from the continued 
fraction representation either directly as poles when the fraction terminates (the high- 
and low-density phases), or using Worpitzsky's theorem (the maximal current phase). 

For a chain of length N the one- and two-point density correlation functions for the 
PASEP may be obtained in terms of the matrices D q and C q as 

(n) = ^(w\Cl- l D q c^\v) , 

(TiTj) = ^{W\Cl- l D q Cp~ l D q C^\V) . (6.64) 

We have already evaluated the thermodynamic limit N — > 00 of by picking up the 
leading singularity in the partial fraction representation. If we switch back to the transfer 
matrix picture of C q , the leading singularity of the generating function z cr is just the 
inverse of the largest eigenvalue of C q , \ cr = l/z cr , which dominates the normalization 

as 

(w\c?\v) = 1 £{w\c?\\)(\\v) 

x 

A 
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~ const. . (6.65) 
For the two-point function, there is a similar expansion [21] 

(w\c i - 1 D q q- i - 1 D q c^\v) 

-X^ 2 (W\X cr }(X cr \D q \X cr }(X cr \D q \X cr }(X cr \V) 
+ \Z~ 2 (W\\ sub )(\ sub \D q \\ cr )(\ cr \D q \\ cr )(\ cr \V) (6.66) 



+ \l~ 2 (j^j (W\X cr )(X cr \D q \X sub )(X sub \D q \X cr )(X cr \V) 
fx \ N ~ j 

+ A^- 2 i-f^j (W\X cr )(X cr \D q \X cr )(X cr \D q \X sub )(X sub \V) , 

where the second largest eigenvalue X sub is the inverse of the sub-dominant singularity 
z sub in the generating function Z. From this it is clear that the ratio X sub /X cr , or 
alternatively z cr /z sub , determines the correlation length. 

If we consider first the high- and low-density phases the sub-dominant singularity is 
separated from the dominant one and both can be read off directly from the continued 
fraction representation of Z. In this case we descend two levels in the continued fraction 
in order to extract both the dominant and sub-dominant singularities, which means 
that we impose c 2 = 0, i.e., $ = l/(aq) or a = l/((3q) depending on the values of the 
parameters. The generating function then truncates to 

Z 2 {a,(3,q,z) = 1 — -z^— 

1 — d z — 



1 — d\z 
1 — d\Z 



- - - , (6.67) 

1 - (d + di)z + (d d 1 - cifz 2 

which allows us to extract both the dominant singularity z cr and subdominant 
singularity z sub for various ranges of the parameters a, (3, q. As we have already noted 
in extracting the normalization from the continued fraction, descending to deeper levels 
still picks up z cr and z sub as the dominant and sub-dominant pole at every deeper level, 
in addition to the further sub-leading singularities. 

This in turn allows us to classify various sub-phases within the high- and low- 
density phases which have different behaviours of the correlation length. For example, 
the low-density phase may be divided into the three sub-phases found in [21]: 
(Ai) aq > 1, aq > (3: 

_ 1 - g _ 1 ~ Q (R 

Zcr (1 + a)(l + a" 1 ) ' Zsub (1 + aq)(l + (aq)-^) ' 

(A 2 ) a > (3 > aq, (3 > 1: 

Zcr = (1 + aXl + a- 1 ) ' ^ = (1 + /)(1 + ^) 5 (6 ' 69) 
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1/P 



A, 



A 3 




Figure 4. The various sub-phases Ai, A2, . . . of the PASEP defined by the behaviour 
of the correlation length are delineated. 



(A 3 ) a > 1 > aq, (3 < 1: 



Zsub = —r- ■ (6.70) 



The 2 slt 6 values which play a role in the A\ and A 2 sub-phases are the (sub- dominant) 
poles of the truncated Z 2 (a, (3,q, z), whereas in the A 3 sub-phase it is the singularity 
determined using Worpitzsky's theorem for the full, untruncated Z. The high-density 
phase splits in a similar fashion with a and (3 interchanged in the various expressions for 
z cr and z su b above. In all these sub-phases the correlation functions decay exponentially 
and the correlation length is given by the log of the ratio of the poles £ = (In (z sub / ' z^))^ 1 . 
The phase structure deduced in this manner by considering the singularities of Z agrees, 
as it should, with that obtained from considering the integral representation of in [21] 
and is shown in figure |4j 



7. (Motzkin) Path Transitions 

In our earlier work on Dyck path representations of the ASEP normalization it was found 
that the high- and low-density phases of the ASEP were reflected in bound phases for the 
lattice paths, whereas the maximal current phase corresponded to an unbound phase. 
The second-order transition lines of the ASEP corresponded to unbinding transition lines 
for one of the two lattice paths which described the ASEP normalization and the first- 
order transition between the high- and low-density phases of the ASEP corresponded 
to a cooperative transition involving both paths. 

The situation for the q < 1 phase structure of the PASEP is analogous. In the 
region of the phase diagram corresponding to the high- and low-density phases of the 
PASEP predominantly one type of horizontal step is bound to the x-axis, reflecting the 
majority of particles or holes in the high- and low-density phases of the PASEP. In the 
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maximal current phase the entropy dominates and paths are unbound. For instance, in 
the low-density A phases, indicated in figure [4j the horizontal on-axis steps weighted 
by 1 + a will make a greater contribution to the partition function than those weighted 
by 1 + /3 and so will dominate the sum. The reverse is true in the high-density phase 
where the steps weighted by 1 + /3 dominate. 

The non-trivial terms in the weights of steps at height n are proportional to q n 
so these steps do not contribute to the energy as N — > oo when q < 1, but only 
entropically. It is therefore not surprising that for q < 1 the PASEP phase diagram is 
essentially identical to the ASEP, with only the position of the transition lines changing 
as q is varied. It is interesting to inquire about the nature of the paths contributing to 
the partition function for q > 1. The expression for Zn in this reverse bias phase was 
calculated in [4], 

> N 

1 a(5 



<N 



A{a,f3; q) (g" 1 ^, 1/5/3; q- 1 ) 



where 



and 



A(a,P;q) 



M(q) 



71 

Inq 



cxp 



q 



(ln/3/5) 2 
4 In q 



q 4 



M{q) + 



1 1 



^ k q k - 1 ' 

k=i H 



(7.71) 



(7.72) 



(7.73) 



The most interesting feature of (7.71) is that one now sees an area-like behaviour 



TV 2 



emerging, whereby ~ q^ 

To see how such two-dimensional behaviour emerges from a model of one- 
dimensional paths when q > 1 consider a tent-shaped path contributing to Z^ with 
N/2 upward steps followed by N/2 downward steps. This has a weight 

(a(3; q) N /2(q;q)N/2 



(q-l) N 

which has a leading term of the form 
(aft- 1) (&P) ' 



q 4 



(q 



(7.74) 



(7.75) 



The reverse bias expression for Zn in (7.71) displays exactly this behaviour. With a 
slight rewriting of the second factor, 

v N 

5/3- 1 \ I \I&P 



Z N ~ A(a,j3;q)(q l af3,l/qaf3;q x ) 




, (7.76) 



which matches the behaviour of (7.75) coming from a "tent" path. Heuristically this 



makes sense since paths which gain in altitude as quickly as possible will donate the 
highest powers of q to Z N and dominate when q > 1. 
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This change in behaviour at q > 1 to an area-like scaling is thus interpreted in 
the lattice path picture as an inflation transition, where the fluctuating path becomes 
the boundary of an inflated vesicle (with particular edge weights favouring a speedy 
ascent) pinned at its ends to the horizontal axis. The matrices D q and E q initially 
denoted the presence or absence of particles on a particular site in the original matrix 
product ansatz. We can associate an up-step and one colour of horizontal step in a path 
with, say, the presence of a particle in the original PASEP steady state. Looking at the 
form of the leading contribution to Zn when q > 1, namely N/2 up-steps followed by 
N/2 down-steps, we can see that this corresponds to a half-filled lattice, which has also 
been observed in simulations of the reverse bias phase and heuristically derived using 
particle/hole symmetry arguments [4]. 



The form of equation (7.71) is reminiscent of the partition function for inflated 
lattice vesicles with a fixed perimeter discussed in [22] . If one models a two-dimensional 
vesicle with a convex polygon of perimeter 2N on a square lattice then the partition 
function is given by 

W-nSl^. ("7) 

where each component square making up the polygon has a weight q and p < 1. The 
summation arises from rectangles and may be extended to run from ±oo with an error 

TV" 2 

of 0(q~~) and the prefactor can be thought of as arising from nibbling off the corners 
of the rectangles with Ferrers diagrams as described in [22]. In the case of the PASEP 
the expression for Zn when q > 1 is suggestive of a similar interpretation; the leading 
"tent" term is decorated with nibbled edges from the fluctuations. 

8. The Case q = 1 

The expression for Zn may also be evaluated in the case of the symmetric simple 
exclusion process (SSEP) when q — 1. If we introduce some more standard notation, 
this time for g-numbers, 



1 — a 

n} q = —^- = l + q + ...q n - 1 , (8.78) 
l-q 



(with [n] q = for n < 0) we can rewrite the Motzkin-path weights as 
J„ = 2[n]l+ (I 

-c. = [»].([»-!].+ g + I) 9 --<m^) , (8.79) 

or, splitting the horizontal weights into two types of steps again in order to make the 
correspondence with bicoloured Motzkin paths explicit, 

q n 



d n ,l = [n)q + 



a 



d„,2=[n] q + j. (8.80) 
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The q- integers become normal integers when q — 1, so the weights at q — 1 are given by 

1 

«n,i =n+- , 
a 

d~n,2 = n + - , 

n(n-l + - + -l , (8.81] 



a /3 

which simplify even further when a = /3=lton + l and n(n + 1). This gives 

2(1,1,1,*) = l — 7 r- 2 , (8.82) 

1 - 2z = 

2 x 3z 2 

1-Az 



3 x Az 2 
l-6z 



which can be seen to be a continued fraction expansion for the divergent power series 

oo 

Z=^2{N + l)\z N , (8.83) 

Af=0 

so Z N — (N + 1)! in this case. A similar expansion exists for general a and (3 and gives 
Zn as the ratio of two Gamma functions, 

7 _ r (A + iV + l) , , 

Zn ~ r(A + i) (8 ' 84) 

where A = ^ + 4 — 1. This agrees with a direct calculation in [5]. Using Stirling's 
approximation for the Gamma functions, and remembering that A is of 0(1) we see that 
InZjv ~ iVlniV so Z^- ~ e NlnN &t q = I, intermediate between the linear behaviour for 
q < 1 and the area law behaviour when q > 1. 

9. The Case q n+1 = 1 

The continued fraction representation of Z will also terminate when q n+1 = 1. These 
values of q have recently proved of interest for the Bethe Ansatz approach to the open 
XX Z spin chain with non-diagonal boundary conditions [23-25]. The Hamiltonian of 
the XX Z spin chain is equivalent (up to a unitary transformation) to the transfer matrix 
of the PASEP. If we denote the expression for A n (z) when q n+1 = 1 by A n (z) the sum 



in (4.48) truncates to two terms 

z n+1 ( a \ 

Ai{z) = _ g )n+i a n+i \S & P> tir+i ~ ^ ' Q) n+ii^- 19 ; q) n+1 j (9.85) 

where we have used a limiting procedure 

i im (r n+1 -,q)n + i = lim ( _ 1)n+1 (n+1)(n+2)/2 = _ x (g g6) 
^° (q; q)n+i 

with q = exp(e + 27rim/(n + 1)) to handle the indeterminate final term [26]. Noting 
that (ae 10 ,q) n+ i = (1 — <5 n + 1 e «( n + 1 ) e ) when q n+1 = 1 this can be written as 

= (1 - q )n + i { ei(n+1)e + e ' l{n+1)(> - + ^" +1 )J ( 9 - 87 ) 
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and we find that the zeros of A n (z) are given by 

Zk = n J. WZ i S (9 - 88) 

(l + nf)(l + ^) 
where r is the root with smaller argument of 



(a n+1 4- B n+1 ) / (a n+1 + 3 n+1 ) 2 
r n+1 = [ - + \l [ - + / ~ 1 ( 9 - 89 ) 

or 



It is interesting that the overall structure of the roots remain similar to that observed 
when a = l/(g n /3), but one no longer has a dominant pole appearing at first order and 
remaining. 



10. Conclusions 



We have shown that the matrix C q = D q +E q in a particular tridiagonal representation of 
the PASEP algebra can be interpreted as the transfer matrix for weighted Motzkin paths, 
with two colours of horizontal steps. Writing the generating function for these paths 
as a continued fraction allowed a succinct derivation of the thermodynamic limit of the 
normalization of the PASEP in its various phases. In particular it allowed calculations 
without explicitly summing the generating function in closed form, which was possible 
(or at any rate, easy to perform) only in the q = case. Consideration of the sub-leading 
singularities in the continued fraction also allowed the determination of the correlation 
length in the high- and low-density phases. 

A further interesting feature of the continued fraction respresentation of the 
generating function was that it made clear how the finite-dimensional representations of 
the PASEP algebra, valid only along special lines in the phase diagram, related to the 
general, infinite-dimensional solution via truncation. The unusual structure of the poles 
of the convergents of the continued fraction along these special lines was highlighted, 
with lower order poles remaining present as the order of the convergent increased and 
the dominant pole (for the high- and low-density phases) already being present at zeroth 
order. The presence of the maximal current phase, on the other hand, was deduced from 
the continued fraction respresentation by using Worpitzky's theorem. 

The phase transitions of the PASEP can be identified with transitions in the lattice 
path model. When q < 1 there are similar unbinding transitions to those seen in the 
Dyck path model for the ASEP, with bound states corresponding to the high- and low- 
density phases of the PASEP and an unbound phase corresponding to the maximal 
current phase. For the reverse bias phase of the PASEP when q > 1 the lattice path 
model was seen to be in an inflated phase and the leading contribution to the lattice 
path partition function was identified with a half-filled state in the PASEP. 
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Interestingly, bicoloured Motzkin paths play a major role in various combinatorial 
bijections and the ASEP and PASEP have recently been investigated for their 
combinatorial interest [27-31]. One common theme has been the appearance of various 
well-known combinatorial weights as special cases for the (P)ASEP normalization 
which can be related to known results in the enumeration of permutations and other 
combinatorial objects. The approach in [29] in particular is rather similar to that 
espoused here, namely going directly to a matrix representation for D and E and seeking 
a combinatorial interpretation (although the principal focus in [29] was on a different 
representation and permutation tableaux). 

It would be an interesting exercise to relate the phase transitions discussed here 
for Motzkin paths and the related phase structure of the PASEP to conformational 
transitions in other (weighted) combinatorial objects such as parallelogram polyominoes, 
binary trees and the permutation tableau of [29] . It would also be worthwhile to relate 
the lattice path picture to discussions of the large deviation functional of the (P)ASEP, 
such as that in [32], where the appearance of a combination of a Brownian excursion 
and a Brownian walk in the discussion is strongly suggestive of a polyomino bijection 
from the Motzkin paths here. 
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Appendix A. Parameters Appearing in the Tridiagonal Matrix 
Representation for 7,5^0 

When all four boundary rates, a, 0, 7 and 5 are nonzero, the tridiagonal D and E 
matrices (5.55) contain a number of parameters. These are 

n-l 



di 



dl 



(1 - q 2n - 2 abcd)(l - q 2n abcd) 

x [bd(a + c) + (b + d)q - abcd(b + d)q n ' 1 - {bd(a + c) + abcd(b + d)}q r ' 
- bd{a + c)q n+1 + ab 2 cd 2 {a + c)q 2n - 1 + abcdib + d)q 2n ] , 



(1 - q 2n ~ 2 abcd)(l - q 2n abcd) 

x [ac(b + d) + (a + c)q — abcd(a + c)q n ~ 1 — {ac(b + d) + abcd(a + c)}q r ' 
- ac{b + d)q n+1 + a 2 bc 2 d{b + d)q 2n ~ 1 + abcd{a + c)q 2n ] , 



1 A 

1 — q n ac 

q n bd 



1 - q n bd 



A 



q ac 
1 — q n ac 
1 

1 - q n bd 



-A 



A 
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which involve the further parameters 

1 



d 



A 



2a I 
1 

2f3 
1 

2a 
1 

2(3 



1 — g — a + 7) + a/(1 - q - a + •y) 2 + 4cry 



[l-q-/3 + S) + v / (l-q-p + 5) 2 + 4:f36 



[1 - q - a + 7) - a/(1 - q - a + y) 2 + 4c*7 
1 - g - /3 + 5) - ^(l-q- (3 + 5) 2 + 4/55 



(l-g n ~ 1 afeed)(l-g"+ 1 )(l-q"afc)(l-g n ac)(l-g' 1 ad)(l-g"fec)(l-i?"fed)(l- 1 ?"cd) 
( 1 - g 2 ™ - 1 a&cd) ( 1 - g 2 ™ afccd) 2 ( 1 - g 2 ™+ 1 abed) 



Finally, the constant appearing in (|5.56|) is 

h 



(abed; q) c 



(q, ah, ac, ad, be, bd, cd; q) c 



1/2 
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